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Abstract 


Twenty-five typical massive white dwarfs (WDs) are selected and the proton decay reaction catalyzed by magnetic 
monopoles (MMs) for these WDs is discussed. A velocity-dependent correction factor strongly affects the cross-section. 
We find that a strong suppression controls the monopole catalysis of nucleon decay by the correction factor. The 
maximum number of MMs is captured and the luminosity can be 2.235 x 10?! and 1.7859 x 10*erg s^! (e.g., for the 
O--Ne core mass WD J055631.17+130639.78). The luminosities of most massive WDs agree well with the 
observations at relatively low temperatures (e.g., 76 = 0.1), but can be three and two orders of magnitude higher than 
those of the observations for model (T) and (ID at relatively high temperatures (e.g., Tg — 10), respectively. The 
luminosities of model (I) are about one order of magnitude higher than those of model (II). Since we consider the effect 
of the number of MMs captured on the mass-radius relation and the suppression of the proton decay by the correction 
factor, the study by model (IT) may be an improved estimation. 
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1. Introduction 


White dwarfs (WDs) form at very high temperatures. Their 
radiation decreases over time from the initial high color 
temperature to turn red. Some observations show that the 
effective temperature of WDs may be mostly between 
5500 K - 40000 K (e.g., Gao et al. 2023), and a few are outside 
this range, and the total thermal energy may be less than 10^ 
erg. What are the specific heat sources in WDs? What can 
provide this energy to WDs? These issues will be discussed in 
this paper. 

From the Stefan-Boltzmann law, the radiation energy 
luminosity is defined by the surface effective temperature 
Тер» 


Lag = АпК2оТА, (1) 


where Lraa and А are the radiation luminosity and radius of the 
WD respectively. с = 5.6704 x 10 ^ ergs ! cm ?K * is the 
radiation constant. 

Mestel (1952), Avakian (1972) and Bildsten & Hall (2001) 
investigated the problem of WD energy sources. Their results 
showed that small amounts of Ne in WDs may be an extra 
source of heat. A possible heat source may be considered from 
22е sedimentation (e.g., Deloye & Bildsten 2002; Garcla- 
Berro et al. 2008). Lobato et al. (2018) and Cheng et al. (2019) 
also discuss the cooling anomaly of high-mass WDs. They 
pointed out that this problem is from ?Ne. By using molecular 
dynamics methods and phase diagrams, Caplan et al. (2020) 


stated that Ne cannot be a possible cause of the heating 
in WDs. 

In this paper, we present a new method to solve the energy 
source problem for 25 typical massive WDs selected from Kilic 
et al. (2021). Our model is based on the magnetic monopole 
(hereafter MM) catalytic proton decay (RC effect) (e.g., 
Rubakov 1981; Callan 1983). Issues related to MMs are 
forefront subjects in astrophysics (e.g., Detrixhe et al. 2011; 
Rajantie 2012; Fujii & Pierre 2015; Frank et al. 2019; 
Kain 2019). We studied the problem of MMs and other related 
issues (e.g., Peng et al. 2017, 2020; Liu 2013, 2014; Liu & 
Gu 2016a; Liu & Liu 2017c, 2018a; Liu 2016b; Liu & 
Liu 2016c; Liu 2016d; Liu et al. 2017a, 2017b; Liu & 
Liu 2018b). Some scholars also discussed other energy 
mechanisms of WDs, such as magnetic field evolution, the 
state equation of the WDs and related issues (e.g., Deng et al. 
2020, 2021; Li et al. 2016; Li & Gao 2023; Li et al. 2023; Fu 
et al. 2020; Gao et al. 2016, 2017, 2019, 2023; Dong et al. 
2023; Wang et al. 2020; Yan et al. 2021; Zhu et al. 2016). Their 
researches provide a database for further exploration of the 
evolution of WDs. Recently, the energy source problem for 
WDs has been discussed by Peng & Liu (2023); Liu et al. 
(2023a, 2023b); Liu & Liu (2023c) and Liu & Liu 2023d). 
However, the influence of the velocity-dependent correction 
factor on the cross-sections in the Rubakov process is ignored 
in the above works. Because the velocity-dependent correction 
factor strongly affects the proton decay thermonuclear reac- 
tions, we re-examine this issue in this paper. 
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The plan of this paper is given as follows. In the next section, 
we present our recent works on the number of possible MMs, 
and the luminosity by MM catalytic proton decay. The MM 
model and the RC luminosity in WDs are shown in Section 3. 
Results and discussions are given in Section 4. Section 5 
summarizes some conclusions. 


2. Our Recent Works on the Numbers of MMs 
Captured and the Luminosity by Catalytic Proton 
Decay in Stars 


Recently, we discussed the energy source problem for WDs. 
First, we discussed the numbers of MMs captured in space at 
the surface of stars, which is given by (Peng et al. 1985; Peng 
& Liu 2023; Liu et al. 2023a, 2023b; Liu & Liu 2023c, 2023d) 


0 

Natu) = 5.7 X 10 55v nd o, (2) 
where v,, is the velocity of the MMs in space, у з = v,,/ 10 2c, 
mo = m,/10^m,, 105m, < mm < 10'5m,, c is the speed of light 
and P is the amount of MMs in space. The superscript (0) 
denotes space. nf = (104—109) cm-? is the number density of 
baryons in the space of the Milky Way (Peng et al. 1985) and the 
Newton saturation value ¢, is the maximum number of MMs, 
which is defined as С, = Gmgm,, fs x 1.9 x 107227" 


m 10?m, ` 
£n = 3hc/2e and mg, which are the magnetic charge of MMs 
and the baryon number mass, respectively. The MM mass is about 
y, ^ 10 5m, (Polyakov 1974), where m, is the mass of a proton. 

The total number of MMs trapped in space after the 
formation of stars (or planets) is estimated as 


Nn (tot) = AnR?IN,, (sur) 


0 DA 
= 3.463 x Wonton SE (E) 


5 


= 3.463 х 107 anf to R2, (3) 


where Ry, = R/Ro, to = 1/10 yr and £= moF/¢,. R, t and F are 
the radius, age of the star and MM flux in space, respectively. 
The reactions рМ — е*т°М + debris(85%), апа pM — 
e' L^ M + debris(1596) are named proton decay reactions cata- 
lyzed by MMs, which were proposed by the RC effect. The 
luminosity due to the RC effect is (Peng et al. 1985; Peng & 
Liu 2023; Liu et al. 2023a, 2023b; Liu & Liu 2023c, 2023d) 


AmNg(OnVt) mgc? = №пв(оһут) твс2, (4) 


Lin 25 А 

z 
where re, п„ and ng are the radius of the stellar central region, 
and the number density of MMs and nucleons, respectively. 
The reaction cross-section of the RC effect is about 
0,7210 75 10 28 ст2. 1mgc? + 1GeV x 1.6 x 107° erg s. 
ут is the thermal movement speed of the nucleus relative to 
the MM, and ут = JKT/mg ~ 8.691 x 10°T'/? cm/s, where 
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T is the temperature, k is the Boltzmann constant and 
mg ~ 1.78 х 10724 is the nucleon mass. The number density 
of nucleons can be written as 


-3 
ng = 2.90 x (2 (Mig) ? = 2.242 x 10^ M,R? ст, 
8 


(5) 


where А, = 2.96 х 105 M, is the Schwarzschild radius, and 
Mi; = M,/10?, M, = M/M., and R, = К/К (Peng et al. 
1985; Peng & Liu 2023; Liu et al. 2023a, 2023b; Liu & 
Liu 2023c, 20234). 

As a general rule, the values of reaction cross-section are 
distributed between 10 ?^ and 10 ?* cm?. Ma & Tang (1983) 
gave an estimation of c,, zz 4.28676 x 10 7^ cm? by using the 
SU(5) Grand Unified Theory (GUT). 


3. The MM Model and RC Luminosity in WDs 
3.1. MM Catalytic Proton Decay Model (I) in WDs 


The mass-radius relation of WDs is one of the interesting 
issues for astrophysicists. For zero-temperature stars, the 
equations of hydrostatic equilibrium and mass conservation 
are given by (Cox & Giuli 1968) 


dP Gm,p 

ae = um X (6) 
aM = Anr?p, (7) 
dr 


where r is the radial variable (when r = 0 at the center) and m, 
is the mass inside a sphere of radius r. 

According to Equations (6)-(7), the mass-radius relation can 
be expressed as the following 


Po 
Ry & ae 375 
oG (Poke) 


МУ = 1.080 x 1075701. 


(8) 


where Ро = «m c^ /(3h°) ~ 5.9637 x 10? dyn/cm?, pe —A/Z 
is the molecular weight per electron, ро = nou zz 9.6838 x 
10° /ст?, по = 8лт?с? /(3h*) ~ 5.8336 x 107? cm ? and и is 
an atomic mass unit. According to Equations (5) and (8) we 
have 


ng = 2242 x 10 M,(R2)? = 1.7794 x 10992 M? om. 
(9) 


According to Equations (3) and (9), in the case without 
considering the RC effect on the mass-radius relation, we can 
estimate the number of MMs captured from space in the 
lifetime of the WD's progenitor star as 


№ = Nagog (D = 3.463 х 107£v anf to (RL). (10) 


According to Equations (4) and (8)-(10) in the case without 
considering the RC effect on the relation of mass-radius, the 
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total luminosity due to the nuclear decay reaction catalyzed by 
MMs is given by 


Ат 3 ж ж 
І„(@ ғ wt ny ng (05, VT) твс2 = Mng (nvr) mgc? 


= 5.541 x 10^n£v and (os v1) Eto(Re). (11) 


3.2. MM Catalytic Nuclear Decay Model (II) in WDs 


When only the RC effect is considered as an energy release 
process, the equation of state in WDs can be obtained by Izawa 
(1986) 


OP | GM (12) 
OM (r) 4nr^ ' 
о, En (13) 
OM (r) 
aM = Arr’, (14) 
Or 


where Em is the energy generation rate per unit mass resulting 
from mass annihilation. 

According to Equations (12)-(14), the mass-radius relation 
is determined by Izawa (1986) 


RP x2 x (a 


oO 


—2.35 
| NU - 2 x їои 735255, 


(15) 


From the number of MMs captured and the RC effect for 
some low-mass WDs, the number of MMs captured in WDs is 
given as 


= МУ уу (ID = 3.463 x 107 anQt((RI PY. (16) 
According to Equations (5) and (15) we have 


MR — 2242 x 107 M, (RMP) 
= 2.803 х 1055 135 M805 cm. (17) 


According to Equations (4) and (15)-(17), in the case that 
considers the RC effect and the correction factor effect, the 
cross-section of the RC effect is corrected, and the total 
luminosity is expressed as 


Lan N Are Cmn (sv) mgc? 
= ММ пк (олут) твс? 
= 5.541 х 10“пМКу Ст? (o, vq) ég(RMRY.. (18) 


In Equation (18), С is the factor describing the correction 
effect, which may correct the cross-section of the RC effect. 
Due to an extra angular momentum carried by the monopole 
electric charge system, Arafune & Fukugita (1983) discussed a 
velocity-dependent correction factor, which will strongly affect 
the cross-section because of the RC effect. The correction 
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factor is given by 


Bo 


where B= v,,/c, Bo =1/(romyA*”?), ro = 1.2 fm and my and А 
are the neutron mass and nuclear mass number, respectively. In 


бе 5 = (Bromy At)’, (19) 


1/2 
the monopole-electric charge system, v = -j + (i + lq D. f 


q=Zx ( > 1,....). 


4. Results and Discussions 


Study of the MM flux has gained considerable interest issue 
in the astrophysical fields since GUTs predicted baryon decay 
reaction and the existence of the MMs (G’t 1974; Parker 1970). 
Parker (1970) gave the MM flux as F< 1075 em?s tsr! 
(Parker 1970). A limit on the flux and the catalysis cross- 
section in neutron stars may be F(ov)_23 < 10721ст 2515г! 
(Kolb & Turner 1984). Freese (1984) estimated the ММ flux іп 
WDs, which may be F(ov). »s < 10 5 em 257! sr !. Freese & 
Krasteva (1999) also showed that the bound was stated as 
F(av). og < 107 28 0m 2515г !. In this paper, we select some 
typical parameters as follows: the number density of baryons 
ng = 10? cm? and the temperature Te = 0.01, 0.1, 1, 10. The 
mass of MM is selected as ту = 10 GeV, and the MM flux is 
F=19x 10 ? cm ?s !sr |. According to the relation 
€=moF/¢,, we have € = 10°. 

Figure 1 shows the correction factor as a function of the 
mass of MM for some typical nuclei due to the RC effect. One 
finds that as the mass of MM increases, the correction factor 
decreases. For example, when т = 90 GeV, the correction 
factors of ^He, эс, 160, 2861, 40 Ar, ^'Ca and ?ÓFe сап be 
calculated as 6.902 x 107, 4.022 x 1077, 2.803 x 10^, 
2.699 x 107^, 8.207 x 107*, 5.239 x 1074 and 2.519 x 107°, 
respectively. These results show that a strong suppression 
controls the monopole catalysis of nucleon decay. Due to the 
effect of these correction factors, the cross-section for the 
Rubakov process also receives a strong suppression. This 
suppression is caused by a long range force, which originates 
from an extra angular momentum carried T a monopole- 
electric-charge system (1.е., q = Z x ( P ...) with Z the 
charge of the nucleus, e.g., Arafune & e 1983). For 
example, for some heavy elements more important suppression 
arises from a strong repulsive force against a slowly moving 
MM because of the Zeeman effect and the diamagnetic effect 
on the atomic electrons (e.g., Malkus 1951). Even in the case of 
light atoms of helium, the repulsive potential is very small (e.g., 
it is about 16 keV for helium). 

Figure 2 shows the number of MMs captured from space in 
the lifetime of the O--Ne core (panels (a)-(d)) and C4-O core 
((e)~(f)) high-mass WDs (Kilic et al. 2021) for model (I) and 
(II) as a function of M,, due to the RC effect when Є = 10? and 
np = 106 cm? at the temperature of Tg = 0.01, 0.1, 1, 10, 
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Figure 2. The number of MMs captured from space in the lifetime of the O+Ne core high-mass WDs ((a)-(d)) and the C4-O core high-mass WDs ((e)-(h)) (Kilic 
et al. 2021) for model (I) and (II) as a function of M, when © = 10? and ng = 106 cm? at the temperature of 76 = 0.01, 0.1, 1, 10, respectively. 


respectively. As the temperature increases, the number of MMs increases by about two orders of magnitude. (ie., it first 
captured increases by about two and three orders of magnitude increases from 3.410 x 107 to 1.078 x 10/5, then increases 
for model (I) and model (ID, respectively. For example, for from 3.410 x 10!8 to 1.078 x 10?) for WD J055631.17 
model (I) in panels (a)-(d), the number of MMs captured +130639.78 (M. = 1.207, T.g= 8340 К, to = 3.33) when the 
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Figure 3. The luminosity as a function of Му for the O--Ne core high-mass WDs ((a)-(d)) and the C4-O core high-mass WDs ((e)-(h)) (Kilic et al. 2021) when 
mg = 10 GeV, C = 7.037 x 107°, E= 10° and ng = 106 cm™°? at the temperature of Tę = 0.01, 0.1, 1, 10, respectively. 


temperature increases from То = 0.01 to Tg = 10. However, the 
number of MMs captured increases by about three orders of 
magnitude for model (П) (i.e., it increases from 3.322 х 10/5 to 
2.352 x 107). 

When the temperature is certain, one sees that the number of 
MMs captured № > №. The higher the mass is, the larger the 
number of MMs captured becomes from Figure 2. For 
example, the maximum MMs captured are N, = 1.078 х 10? 
and № = 2.352 x 10?!, for the O+Ne core high-mass WD 
J055631.174-130639.78 for model (I) and (II) at Tg— 10, 
respectively. However, the maximum MMs captured are 
№ = 1.053 x 10? and № = 1.907 x 10?!, for C+O core 
high-mass WD J055631.174-130639.78. On the one hand, by 
comparing the number of MMs captured of low-mass C+O 
core WD J055631.17--130639.78 with that of high-mass C+O 
core WD J190132.74+145807.18, one can find that the 
number of MMs captured of WD J055631.174-130639.78 is 
М, = 1.053 x 101°, which is larger than that of WD J190132.74 
+145807.18 (i.e., №, = 1.016 x 10!%) for model (I). Under the 
same condition, we can see from Figure 2 that the numbers of 
MMs captured аге N; = 1.907 x 107! and № = 1.365 x 10'? 
for WD J055631.17+130639.78 and WD 119013274 
+145807.18, respectively for model (П) at Tg = 10. The main 
reason is that according to Equations (8) (10), (15) and (16), we 
have Ni оос 1/M., but Л 2 « Ts. On the other hand, the effect 
of mass-radius on MMs captured is ignored for model (I). 
Model (ID) considers the effect of mass—radius and the RC 
effect on the number of MMs captured, so the data are 
relatively accurate. 


The mass-radius relation in model (I) is considered only on 
the condition of zero temperature according to classical formula 
Equation (8). According to Equation (15), the mass-radius 
relation in model (II) is considered on the condition of zero 
temperature, as well as the RC effect and the number of MMs 
captured N, only for some lower mass WDs. It is found that, if 
a star contains a sufficient number of monopoles, its radius 
increases after some stages and finally it dissolves into a diffuse 
state. However, its lifetime is much longer than the age of the 
universe because an energy source is rest mass, which is much 
larger than the nuclear energy. 

From Figure 1, when mg = 10 GeV, the correction factors of 
"C and '°O can be 7.037 х 107, and 8.663 x 10^, 
respectively. However, the correction factors can Бе 
6.821 x 107^ and 5.279 x 107^ when т = 60 GeV. In this 
paper, we select т = 10 GeV and the correction factor 
C — 7.037 x 10? for our WD samples (Kilic et al. 2021). 
Figure 3 displays the luminosities as a function of Му for WDs 
(Kilic et al. 2021) at the temperature of Т = 0.01, 0.1, 1, 10 
when mjg—10GeV, C=7.037x10* and ¿=10°, 
nj = l06cm-?. One sees that the luminosities for WDs 
increase as the temperature increases for model (T) and (ID, 
but decrease for model (II). For example, from Figure 3 for 
model (I), luminosities increase from 1.687 x 102° to 
1.687 x 107, for WD J055631.17--130639.78 when the 
temperature increases from Ty = 0.01 to Tg = 10, and increase 
from 5.404 x 1078 to 1.142 x 10?! for model (II). To compare 
the luminosities of model (I) with those of model (ID, one can 
see that the luminosities of model (I) agree well with those of 
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Figure 4. The scaling factor (k; — ES (D/Lyaa; К = EMS D/L; Кз = Lm (CD). VEMR (ID) as a function of M, for the O+Ne core high-mass WDs ((a)-(d)) and the C 


+O core high-mass WDs ((e)-(h)) (Kilic et al. 2021) when т = 10 GeV, C = 7.037 x 1073, ё = 10° апа ng = 106 ста? at the temperature of Tg = 0.01, 0.1, 1, 10, 
respectively. 


Table 1 
Some Information on the High-mass WD Model (Kilic et al. 2021). All solutions above 1.29 Mo for the O+Ne core models and above 1.334 Mo for the C+O core 
models are extrapolated 


O--Ne соге O+Ne core C+0 core C--O core 
Star Star Name Composition Тек Мазз Cooling Age Mass Cooling Age 
(K) (Ms) (Gyr) (Mo) (Gyr) 
1 J010338.56—052251.96 H 9040 + 70 1.262 + 0.003 2.84 + 0.03 1.310 + 0.003 2.60 + 0.04 
2 J025431.454-301935.38 logH/He = — 5 11060 + 560 1.302 + 0.024 2.25 + 0.10 1.330 + 0.016 1.49 + 0.17 
3 J114012.81+232204.7 H 11 860 + 220 1.294 + 0.008 2.10 + 0.04 1.336 + 0.006 1.71 + 0.06 
4 J132926.04+254936.4 H 29 010 + 750 1.314 + 0.006 0.81 + 0.05 1.351 + 0.006 0.37 + 0.03 
5 J172736.28+383116.9 H 9420 + 200 1.252 + 0.012 2.78 + 0.08 1.302 + 0.011 2.59 + 0.12 
6 J183202.83--085636.24 He 34 210 + 1020 1.301 + 0.006 0.45 + 0.03 1.319 + 0.004 0.20 + 0.02 
7 J190132.74--145807.18 H 29 100 + 480 1.279 + 0.003 0.61 + 0.02 1.319 + 0.004 0.35 + 0.02 
8 J221141.80--113604.5 H 9020 + 160 1.262 + 0.009 2.85 + 0.07 1.310 + 0.008 2.61 +0.11 
9 J225513.48--071000.9 H 10 990 + 210 1.252 + 0.012 2.36 + 0.05 1.302 + 0.011 2.18 + 0.09 
10 J235232.30—025309.2 H 10 680 + 100 1.272 + 0.003 2.38 + 0.02 1.319 + 0.003 2.10 + 0.03 
11 J004917.14—252556.81 H 13 020 + 460 1.263 + 0.011 1.94 + 0.08 1.312 + 0.010 1.72 + 0.09 
12 Ј032900.79—212309.24 Н 10 330 + 290 1.305 + 0.010 2.32 + 0.06 1.344 + 0.008 1.87 + 0.09 
13 J042642.02—502555.21 H 17 900 + 1570 1.264 + 0.019 1.30 + 0.16 1.312 + 0.016 1.08 + 0.16 
14 J043952.72+454302.81 H 19 120 + 630 1.258 + 0.008 1.18 + 0.06 1.307 + 0.007 0.96 + 0.06 
15 J055631.17--130639.78 H 8340 + 260 1.207 + 0.021 3.33 + 0.12 1.257 + 0.023 3.34 + 0.18 
16 J060853.60—451533.03 H 19 580 + 1910 1.258 + 0.021 1.13 + 0.16 1.307 + 0.019 0.92 + 0.17 
17 J070753.00--561200.25 H 18 100 + 350 1.240 + 0.005 1.23 + 0.04 1.291 + 0.005 1.06 + 0.04 
18 J080502.93—170216.57 H 10 830 + 110 1.254 + 0.004 2.40 + 0.03 1.304 + 0.003 2.20 + 0.03 
19 J093430.71—762614.48 H 10 050 + 1350 1.284 + 0.055 2.47 + 0.35 1.328 + 0.047 2.11 + 0.50 
20 J095933.33—182824.16 H 12 000 + 180 1.273 + 0.005 2.12 + 0.03 1.320 + 0.004 1.83 + 0.04 
21 J111646.44—160329.42 H 10 480 + 170 1.264 + 0.007 2.45 + 0.05 1.312 + 0.006 2.21 + 0.07 
22 J125428.86—045227.48 H 14 420 + 390 1.258 + 0.008 1.71 + 0.06 1.308 + 0.007 1.52 + 0.06 
23 J174441.56—203549.05 H 27 140 + 890 1.271 + 0.008 0.65 + 0.06 1.312 + 0.008 0.43 + 0.04 
24 J180001.21+451724.7 H 16 410 + 290 1.253 + 0.003 1.44 + 0.03 1.303 + 0.004 1.26 + 0.04 
25 J181913.36—120856.44 H 37 970 + 1940 1.305 + 0.009 0.37 + 0.03 1.327 + 0.006 0.12 + 0.03 
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Table 2 
The Comparisons Among Several Luminosities Due to the RC Effect for the O--Ne Core High Mass WD model (Kilic et al. 2021) at Relatively Low Temperature 


Ts = 0.1, in which the Parameters 1.207 < Му < 1.315, тє = 10 GeV, C = 7.037 x 102, € = 10? and ng = 106 cm™ (note ky = Lin(D)/Lraas k2 = LMR OD ИТ. аа, 
and ks = Lm (1)/Lm'* (ID) 


Star M. Тек to М, № Г.а Lin) LMR ap kı ky Кз 

1 1.262 9040 2.84 8.9280е16 1.7278е19 2.2800е29 1.5266e30 3.4703e29 6.6926 1.5214 4.399] 
2 1.302 11060 2.25 6.9280e16 9.4080e18 5.0100e29 1.2608e30 3.4214e29 2.5189 0.6835 3.6852 
3 1.294 11860 2.10 6.4920e16 8.5206e18 6.6500e29 1.167230 3.2207e29 1.7562 0.4846 3.6239 
4 1.314 29010 0.81 2.4790e16 1.2860e18 2.3548e31 4.5950e29 2.0170e29 0.0195 0.0086 2.2780 
5 1.252 9420 2.78 8.7860e16 1.7305e19 2.7000e29 1.4786e30 3.3316e29 5.4687 1.2322 4.4381 
6 1.310 34210 0.45 1.3860e16 4.4400e17 4.5842e31 2.5190e29 1.4078e29 0.0055 0.0031 1.7893 
7 1.279 29100 0.61 1.9010e16 8.6600e17 2.4275e31 3.3380e29 1.5638e29 0.0138 0.0064 2.1345 
8 1.262 9020 2.85 8.9600e16 1.7393e19 2.2600e29 1.5320e30 3.4770e29 6.7759 1.5379 4.4060 
9 1.264 10990 2.36 7.4110e16 1.2052e19 4.9800e29 1.2713e30 3.1525e29 2.5542 0.6334 4.0326 
10 1.272 10680 2.38 7.4430e16 1.1844e19 4.4200e29 1.2929e30 3.2409e29 2.9248 0.7332 3.9892 
11 1.263 13020 1.94 6.0960e16 8.3400e18 9.8100e29 1.0439e30 2.822229 1.0641 0.2877 3.6990 
12 1.315 10330 2.32 7.1320e16 9.8510e18 3.8000e29 1.3041e30 3.5089e29 3.4288 0.9226 3.7165 
13 1.266 17900 1.30 4.0780e16 3.8490e18 3.4990e30 7.0180е29 2.2842е29 0.2006 0.0653 3.0723 
14 1.258 19120 1.18 3.7170е16 3.3120е18 4.5740е30 6.3160е29 2.1161е29 0.1381 0.0463 2.9848 
15 1.207 8340 3.33 1.0784е17 2.9607е19 1.7000е29 1.6868е30 3.2192е29 9.9089 1.8911 5.2396 
16 1.279 19580 1.13 3.5210е16 2.7940е18 4.9750е30 6.1840е29 2.1951е29 0.1243 0.0441 2.8170 
17 1.240 18100 1.23 3.9120е16 3.8680е18 3.7090e30 6.4580e29 2.0541e29 0.1741 0.0554 3.1443 
18 1.254 10830 2.40 7.5770e16 1.2978e19 4.720029 1.2792e30 3.0908e29 2.7110 0.6550 4.1387 
19 1.284 10050 2.47 7.6760e16 1.2093e19 3.4400e29 1.3587e30 3.4231e29 3.9445 0.9938 3.9691 
20 1.278 12000 2.12 6.6090e16 9.2730e18 7.0200e29 1.1589e30 3.0938e29 1.6501 0.4405 3.7458 
21 1.271 10480 2.45 7.6660e16 1.2567e19 4.1000e29 1.3295e30 3.2836e29 3.2422 0.8008 4.0490 
22 1.266 14420 1.71 5.3640e16 6.4800e18 1.4740e30 9.2310e29 2.6559e29 0.6264 0.1802 3.4756 
23 1.275 27140 0.65 2.0290e16 9.9300e17 1.8405e31 3.5420e29 1.6010e29 0.0192 0.0087 2.2124 
24 1.253 16410 1.44 4.5490e16 4.9380e18 2.4890e30 7.6670е29 2.3270е29 0.3081 0.0935 3.2948 
25 1.305 37970 0.37 1.1370е16 3.0100е17 6.9426е31 2.0800е30 1.2784е29 0.0030 0.0018 1.6269 


From the value of Кз, the luminosities for model (I) can be only 
3.121 times larger than those of model (П), but can Бе 14.77 
times that at 7g = 10 for ће C+O core WD J055631.17 
+130639.78 from Figure 4(a)-(d). The same conclusions may 
be obtained from Figure 4(e)-(h). 

Tables 1—3 show the comparisons of the luminosities and the 
scaling factor k due to the RC effect for these high-mass O+Ne 
соге WDs at Tg = 0.1, 10, respectively. The luminosities for 
most of the O4-Ne core WDs for model (I) and model (II) agree 
well with the observations and the difference among them is no 
more than one order magnitude at 7 = 0.1. However, as 
temperature increases (e.g., Т = 10 in Table 3), the observa- 
tions can be three and two orders magnitude lower than those 
of model (I) and model (ID, respectively. The scaling factors ky, 
kə can be in the range of 0.0030<k, < 9.9089 and 
0.0018 < k2 < 1.8911 for the O-+Ne core star WD 
J055631.174-130639.78 and WD J181913.36-120856.44, 
respectively at 75 = 0.1 in Table 2, but are in the range of 


model (II) at the relative temperature (e.g., Te = 0.01), but as 
the temperature increases, the difference between the two 
models becomes bigger and bigger. To compare the luminos- 
ities of observation L,,q with those of model (I) and model (ID, 
we find that the difference of the luminosities between the 
observation Laqg and our results can be no more than three 
orders magnitude. For example, for the high mass C4-O core 
WD J190132.74+145807.18, the difference is no more than 
one order magnitude at Tg = 0.01, but, for the C+O core WD 
J055631.17--130639.78, the luminosities of the observation 
Lyaa are about two and three orders magnitude lower than those 
of model (П) and (T), respectively. 

Figure 4 displays the scaling factor k (k= 
Lm Q)/Lrags k2 = La ^ (ID/Lua; k3 = Lm(D/Lm (ID) as a 
function of M, for the O--Ne(C4-O) соге WDs (Kilic et al. 
2021) when mig— 10 GeV, C— 7.037 x 107°, E= 10° and 
np = 106 сш? at the temperatures of Tg — 0.01, 0.1, 1, 10. 


One can see that model (I) and model (II) agree well with the 
observations at relatively low temperature (e.g., Tg — 0.01) 
from Figure 4. As the temperature increases, the difference 
becomes bigger and bigger. For example, from Figure 4, А, 
increases from 3.121 to 990.9 for the C+O core WD 
J055631.17--130639.78, but increases from 0.3175 to 67.10. 


0.2996 < kı < 990.8927 and 0.0653 < ky < 67.1005 at Т, = 10 
in Table 3. Comparing the results from model (I) and model 
(D, one finds 1.6269 < k, < 5.2396 and 4.5851 < k3 < 
14.7673 in Tables 2 and 3, respectively. Table 3 lists 14 
WDs, whose values of k3 exceed a factor of 10 at Tg = 10. For 
example, the values of kz are 12.3983, 10.3862, 10.2135, 


Research in Astronomy and Astrophysics, 24:035014 (11pp), 2024 March 


Liu & Liu 


Table 3 
The Comparisons Among Several Luminosities Due to the RC Effect for the O4-Ne Core High-mass WD model (Kilic et al. 2021) at Relatively High Temperature 
Т, = 10, in Which the Parameters 1.207 < M, < 1.315, тц = 10 GeV, C = 7.037 x 1077, € = 10? and ng = 10° стт? (note kı = Lm(I)/Lrad, ko = LYR dD аа, 
and & = Ln(D)/LM® (ID) 


Star Mx Ter to М, № һа) LY ap kı ky Кз 

1 1.262 9040 2.84 8.9280е18 1.3724е21 2.2800е29 1.5266e32 1.2313e31 669.2608 53.9800 12.3983 
2 1.302 11060 2,25 6.9280e18 7.4730e20 5.0100e29 1.2608e32 1.2140e31 251.8930 24.2527 10.3862 
3 1.294 11860 2.10 6.4920e18 6.7730e20 6.6500e29 1.1672e32 1.1428e31 175.6229 17.1952 10.2135 
4 1.314 29010 0.81 2.4790e18 1.0220e20 2.3548e31 4.5950e31 7.1570е30 1.9513 0.3039 6.4204 
Э 1.252 9420 2.78 8.7860е18 1.3746е21 2.7000е29 1.4786e32 1.1821e31 546.8672 43.7209 12.5082 
6 1.310 34210 0.45 1.3806e18 3.5300e19 4.5842e31 2.5190e31 4.9950e30 0.5495 0.1090 5.0430 
7 1.279 29100 0.61 1.9001e18 6.8800e19 2.4275e31 3.3380e31 5.5490e30 1.3751 0.2286 6.0160 
8 1.262 9020 2.85 8.9600е18 1.3816e21 2.2600e29 1.5320e32 1.2337e31 677.5939 54.5658 12.4179 
9 1.264 10990 2.36 7.4110e18 9.5730e20 4.9800e29 1,2713e32 1.1185e31 255.4156 22.4730 11.3654 
10 1.272 10680 2.38 7.4430e18 9.4080e20 4.4200e29 1.2929e32 1.1499e31 292.4814 26.0142 11.2431 
11 1.263 13020 1.94 6.0960e18 6.6250e20 9.8100e29 1.0439e32 1.0014e31 106.4135 10.2074 10.4252 
12 1.315 10330 2.32 7.1320e18 7.8250e20 3.8000e29 1.3041e32 1.2450e31 342.8781 32.7349 10.4744 
13 1.266 17900 1.30 4.0780e18 3.0580e20 3.4990e30 7.0180e31 8.1040e30 20.0554 2.3162 8.6588 
14 1.258 19120 1.18 3.7170e18 2.6310e20 4.5740e30 6.3160e31 7.5080e30 13.8078 1.6414 8.4122 
15 1.207 8340 3.33 1.0784e19 2.3517e21 1.7000e29 1.6868e32 1.1422e31 990.8927 67.1005 14.7673 
16 1.279 19580 1.13 3.5210e18 2.2190e20 4.9750e30 6.1840e31 7.7880е30 12.4280 1.5653 7.9395 
17 1.240 18100 1.23 3.9120е18 3.7200е20 3.7090e30 6.4508e31 7.2880е30 17.4127 1.9649 8.8618 
18 1.254 10830 2.40 7.5770е18 1.0309е21 4.7200е29 1.2792e32 1.0967e31 271.0969 23.2414 11.6644 
19 1.284 10050 2.47 7.6760e18 9.6060e20 3.4400e29 1.3587e32 1.2146e31 394.4526 35.2615 11.1865 
20 1,278 12000 2.12 6.6090e18 7.3660e20 7.0200e29 1.1589e32 1.0977e31 165.0079 15.6300 10.5571 
21 1.271 10480 2.45 7.6660e18 9.9820e20 4.1000e29 1.3295e32 1.1651e31 324.2232 28.4119 11.4115 
22 1.266 14420 1.71 5.3640e18 5.1470e20 1.4740e30 9.2310e31 9.4230e30 62.6375 6.3944 9.7956 
23 1.275 27140 0.65 2.0290e18 7.8900e19 1.8405e31 3.5420e3 1 5.6810e30 1.9245 0.3086 6.2355 
24 1.253 16410 1.44 4.5490e18 3.9202e20 2.4890e30 7.6670e3 1 8.2560e30 30.8079 3.3176 9.2861 
25 1.305 37970 0.37 1.1370e18 2.3900e19 6.9426e31 2.0800e3 1 4.5360e30 0.2996 0.0653 4.5851 


12.5082, 12.4179, 11.3654, 11.2431, 10.4252, 10.4744, 
14.7673, 11.6644, 11.1865, 10.5571 and 11.4115 for the O 
+Ne core high-mass WD No.1-3, 5, 8-12, 15 and 18-21, 
respectively. 

Tables 4 — 5 also show the comparisons of the luminosities 
and the scaling factor k for ће C+O core WDs at Tg = 0.1, 10. 
For most of the С+О core WDs, the results of model (1) and 
model (II) agree well with the observations of Г, „а, and the 
difference among them is no more than one order of magnitude 
at Tg = 0.1. However, as the temperature increases (e.g., 
Tg = 10 in Table 5), the observations of Г, „а can be three and 
two orders of magnitude lower than those of model (T) and 
model (ID), respectively. The scaling factors kı and kz can be in 
the range of 0.0008 < kı < 6.4014 and 0.0009 < kz < 1.6509 
for the C+O core WD J010338.56-052251.96 and WD 
J181913.36-120856.44, respectively at Tg — 0.1 in Table 4, 
but the ranges аге 0.0823 < kı < 640.1393 and 0.0310 < о < 
58.5767 at Tg = 10 in Table 5. From model (I) and model (ID, 
ks is in the range of 0.9429 < k; < 4.7758 and 
2.6575 < Кз < 13.4601, shown in Tables 4 and 5, for the C 
+O core WD J055631.174-130639.78 and WD J181913.36- 
120856.44, respectively. 

According to stellar evolution theory, radiative transfer of 
energy is negligible compared with thermal conduction in the 


highly degenerate core of WDs, in spite of the reduced opacity 
of a degenerate gas. It will be shown that the thermal 
conduction is so great that, during the cooling of a WD, the 
temperature is always very nearly uniform in the core. MMs 
trapped inside WDs can catalyze decay of nucleons and 
provide an additional source of internal heat for the star. An 
MM which passes through a WD can easily lose enough energy 
to be captured. Ahlen & Kinoshita (1982) calculated this 
energy source, which comes from the electronic interactions 
and is given by dE/dx ғ 10008 GeV/cm, where p is the 
density of a WD (in g cm >), and 3=v,,/c is the ratio of the 
velocity of the MM to the speed of light as it passes through the 
WD. Far from the WD, the velocity of a typical MM in the 
Galaxy is the virial velocity (1.е., vm = 107° c). When an MM 
falls into the gravitational potential well of a WD, it is 
accelerated and (3 increases to (v2 + 2GM/R)!/? at the surface 
of the WD. The energy loss in traveling through the WD can be 
about 5 х 10!” MeV. Once the MM is captured, it sinks toward 
the center of the WD. If we consider the motion of an MM as a 
harmonic oscillator with a dE/dX damping term, then the 
timescale for the MM to fall from rest to the center can be 
estimated to be about 1000 s. 

Based on the above analysis, we can conclude that with the 
increase in the number of MMs captured by a WD, the 
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Table 4 
The Comparisons Among Several Luminosities Due to the RC Effect for the C+O Core High-mass WD Model (Kilic et al. 2021) at Relatively Low Temperature 
Tg = 0.1, in which the Parameters 1.257 < М, < 1.422, т = 10 GeV, C = 7.037 x 102, = 10° and ng = 10$ cm? (note 
ky = Lm (D/Lraa, ko = LV ID/Lag and ky = Ly (МЕШ) 


Star M. Тек to М, N5 Lia Lin LMR ap kı ky Кз 

1 1.310 9110 2.60 7.9730е17 1.1987е19 2.2900е29 1.4689е30 3.7883е29 6.4014 1.6509 3.8775 
2 1.330 11620 1.49 4.5230е17 3.8400е18 6.0100е29 8.5900e29 2.9477e29 1.4286 0.4902 2.9141 
3 1.336 12080 1.71 5.1750e17 4.8720e18 7.0000е29 9.9170е29 3.2323е29 1.4164 0.4616 3.0682 
4 1.351 29760 0.37 1.1120е17 2.5100е17 2.5601е31 2.1780е29 1.4507е29 0.0085 0.0057 1.5014 
5 1.302 9420 2.59 7.9750e17 1.2292e19 2.6300e29 1.4514e30 3.6967e29 5.5100 1.4034 3.9261 
6 1.319 34410 0.20 6.1000e16 8.8000e16 4.6495e31 1.1400e29 9.4760e28 0.0025 0.0020 1.2033 
7 1.422 30200 0.35 1.0160e17 1.7200e17 2.6238e31 2.2060e29 1.6963e29 0.0084 0.0065 1.3005 
8 1.314 9180 2.61 7.9870e17 1.1881e19 2.3600e29 1.4806e30 3.8388e29 6.2704 1.6258 3.8568 
9 1,313 11150 2.18 6.6750e17 8.4730e18 5.1400e29 1.2354e30 3.467329 2.4028 0.6744 3.5630 
10 1.329 10780 2.10 6.3780e17 7.4020e18 4.4600e29 1.2094e30 3.5503e29 2.7141 0.7967 3.4066 
11 1.322 13480 1:72 5.2420e17 5.2090e18 1.0930e30 9.8360e29 3.1204e29 0.8996 0.2854 3.1523 
12 1.344 10620 1.87 5.6370e17 5.5950e18 4.1700e29 1.0932e30 3.4701e29 2.6241 0.8329 3.1504 
13 1.328 19470 1.08 3.2820e17 2.1010e18 4.7440e30 6.2140e29 2.4560e29 0.1310 0.0518 2.5300 
14 1.307 19750 0.96 2.9480e17 1.8270e18 5.0770e30 5.4070e29 2.1718e29 0.1065 0.0428 2.4897 
15 1.257 10940 3.34 1.0528e18 2.4012e19 4.9100e29 1.7859e30 3.7395e29 3.6407 0.7623 4.7758 
16 1.326 21490 0.92 2.7980e17 1.5610e18 7.0480e30 5.2830e29 2.2364e29 0.0749 0.0317 2.3621 
17 1.291 18450 1.06 3.2820e17 2.3550e18 3.8980e30 5.8730e29 2.1927e29 0.1507 0.0562 2.6785 
18 1.304 10830 2.20 6.7670e17 8.9420e18 4.6000e29 1.2354e30 3.3984e29 2.6872 0.7392 3.6351 
19 1.328 10160 2.11 6.4120e17 7.4990e18 3.5200e29 1.2140e30 3.5498e29 3.4509 1.0091 3.4199 
20 1.320 13350 1.83 5.5830e17 5.9080e18 1.0530e30 1.0444e30 3.2108e29 0.9920 0.3050 3.2528 
21 1.325 10660 2.21 6.7260e17 8.2870e18 4.2700e29 1.2677e30 3.6114e29 2.9691 0.8459 3.5102 
22 1.318 14590 1.52 4.6420e17 4.1850e18 1.5040e30 8.6580e29 2.8832e29 0.5758 0.1918 3.0027 
23 1.320 28030 0.43 1.3120e17 3.7700e17 2.0462e31 2.4540e29 1.4477e29 0.0120 0.0071 1.6952 
24 1.303 16600 1.26 3.8780e17 3.1140e18 2.5390e30 7.0680е29 2.4942е29 0.2784 0.0982 2.8338 
25 1.327 39910 0.12 3.6500е16 3.2000е16 8.3800е31 6.9000е28 7.3150е28 0.0008 0.0009 0.9429 


luminosity of MMs that catalyze nuclear decay increases 
linearly with time until they become the main contribution to 
the total luminosity. One can even observe that for some of the 
oldest WDs, the luminosity may have passed its minimum, and 
some reheating may have occurred. 

For example, for some low-mass WDs with a low lifetime 
(e.g., the O+Ne core WD J055631.17+130639.78 with 
Мұ = 1.207, to =3.33 and Т. = 8340), the luminosity may 
increase linearly with time as the number of MMs captured 
increases from 1.6868 x 10?? erg to 1.6868 x 10°? erg at 
Т,= 0.1 ~ 10 for model (D, but from 3.2192 x 107 erg to 
1.1422 x 10?! erg for model (ID. For some high-mass С+О 
core WDs with a low lifetime, such as WD J190132.74 
+145807.18 (М, = 1.422, tg = 0.35, Т. = 30200), the lumin- 
osity may increase linearly with time as the number of MMs 
captured increases from 2.2060 х 10° erg to 2.2060 x 
102° erg at Tg — 0.1 ~ 10 for model (I), but from 1.6963 x 
10° erg to 6.0190 х 10?? erg for model (Ш). 

Since we have neglected the effect on the mass-radius 
relation by the number of MMs captured by a WD, and only the 
mass-radius relation and the RC effect are considered in model 
(D, the luminosities can be overestimated, compared with 
model (II) (see Equations (8) and (11)). By considering the 
effect on the mass-radius relation by the number of MMs 


captured by a WD (see Equations (15), (18) and (19)), as well 
as the RC effect and the correction factor suppression, one sees 
that the calculations in model (II) agree well with the 
observations, and may be an improved estimation compared 
to model (I). 

The annihilation of magnetic and antimagnetic monopoles 
can cause a significant reduction in the number of MMs and the 
catalytic luminosity of the monopole in WDs. Dicus et al. 
(1982) calculated the annihilation cross-sections of MMs and 
anti-MMs caused by two-body and three-body recombinations. 
Their results showed that this kind of annihilation has little 
effect on the flux and luminosity of MMs. The magnetic field 
may keep the MM and anti-MM distributions far enough apart 
that annihilation is negligible. 


5. Conclusions 


Proton decay catalyzed by MMs plays a key role in the 
evolution of burning WDs. Equations (15) and (16) show that if 
a WD contains enough monopoles, its radius will increase 
during the evolution. However, it may live much longer than 
the age of the universe, so the question of its energy source has 
been the focus of many astrophysical studies. In this work, we 
have investigated the problem of the energy source applied to 
25 massive WDs. Some of our main results are summarized as 


Research in Astronomy and Astrophysics, 24:035014 (11pp), 2024 March Liu & Liu 


Table 5 
The Comparisons Among Several Luminosities Due to the RC Effect for the C--O Core High-mass WD Model (Kilic et al. 2021) at Relatively High Temperature 
Ts = 10, in which the parameters 1.257 < Mẹ < 1.422, тє = 10 GeV, C = 7.037 x 107°, € = 10° and ng = 106 cm? (note kı = Ly, (D/Lraas ko. = LMR(ID/Lig 
and k; = Ly, (D/L (ID) 


m 


Star Mx Ter to М, № Гаа Ln) LY (ШП) kı ky Кз 

1 1.310 9110 2.60 7.9730e18 9.5220e20 2.2900e29 1.4689e32 1.3442e31 640.1393 58.5767 10.9282 
2 1.330 11620 1.49 4.5230e18 3.0510e20 6.0100e29 8.5900e31 1.0459e31 142.8555 17.3939 8.2130 
3 1.336 12080 1.71 5.1750e18 3.8700e20 7.0000e29 9.9170e31 1.1469e31 141.6359 16.3790 8.6474 
4 1.351 29760 0.37 1.1120e18 1.9900e19 2.5601e31 2.1780e31 0.5147e30 0.8508 0.2011 4.2315 
5 1.302 9420 2.59 7.9750е18 9.7640e20 2.6300e29 1.4514e32 1.3117e31 550.9982 49.7956 11.0652 
6 1.319 34410 0.20 6.1000e17 7.0000e18 4.6495e31 1.1400e31 3.3620e30 0.2453 0.0723 3.3915 
7 1.422 30200 0.35 1.0160e18 1.3700e19 2.6238e31 2.2060e31 6.0190e30 0.8408 0.2294 3.6652 
8 1.314 9180 2.61 7.9870e18 9.4370e20 2.3600e29 1.4806e32 1.3621e31 627.0361 57.6847 10.8701 
9 1.313 11150 2.18 6.6750e18 6.7310e20 5.1400e29 1.2354e32 1.2302e31 240.2805 23.9278 10.0419 
10 1.329 10780 2.10 6.3780e18 5.8792e20 4.4600e29 1.2094e32 1.2597e31 271.4103 28.2686 9.6011 
11 1.322 13480 1.72 5.2420e18 4.1380e20 1.0930e30 9.8360e31 1.1072e31 89.9630 10.1260 8.8843 
12 1.344 10620 1.87 5.6370e18 4.4440e20 4.1700e29 1.0932e32 1.2312e31 262.4062 29.5538 8.8789 
13 1.328 19470 1.08 3.2820e18 1.6690e20 4.7440e30 6.2140e31 8.7140e30 13.0975 1.8368 7.1305 
14 1.307 19750 0.96 2.9480e18 1.4520e20 5.077030 5.4070e31 7.7060e30 10.6509 1.5179 7.0168 
15 1.257 10940 3.34 1.0528e19 1.9074e21 4.9100e29 1.7859e32 1.3268e31 364.0657 27.0477 13.4601 
16 1.326 21490 0.92 2.7980e18 1.2400e20 7.0480е30 5.2830е31 7.9350е30 7.4948 1.1258 6.6574 
17 1.291 18450 1.06 3.2820е18 1.8710е20 3.8980е30 5.8730е31 7.7800е30 15.0662 1.9958 7.5492 
18 1.304 10830 2.20 6.7670е18 7.1030е20 4.6000е29 1.2354е32 1.2058е31 268.7176 26.2286 10.2452 
19 1.328 10160 2,11 6.4120е18 5.9560е20 3.5200е29 1.2140е32 1.2595e31 345.0927 35.8037 9.6385 
20 1.320 13350 1.83 5.5830e18 4.6930e20 1.0530e30 1.0444e32 1.1392e31 99.1987 10.8204 9.1678 
21 1.325 10660 2.21 6.7260е18 6.5830е20 4.2700е29 1.2677e32 1.2814e31 296.9117 30.0120 9.8931 
22 1.318 14590 1.52 4.6420e18 3.3250e20 1.5040e30 8.6580e31 1.0230e31 57.5818 6.8041 8.4628 
23 1.320 28030 0.43 1.3120e18 2.9900e19 2.0462e31 2.4540e31 5.1370e30 1.1994 0.2510 4.7777 
24 1.303 16600 1.26 3.8780e18 2.4730e20 2.5390e30 7.0680e3 1 8.8500e30 27.8393 3.4857 7.9867 


25 1.327 39910 0.12 3.6500e17 2.6000e18 8.3800e31 6.9000e30 2.5950e30 0.0823 0.0310 2.6575 


follows: First, by considering the effect of temperature on the respectively (e.g., WD J055631.17+130639.78). Based on 


cross-section, thermonuclear reaction and mass-radius relation, the above analysis and discussion, the monopole-catalyzed 
we presented two MMs models (I, II) to discuss the issue of nucleon decay process could be preventing WDs from cooling 
energy source. Second, we evaluated the number of MMs down into a stellar graveyard by keeping them hot. 

captured by WDs and compared the luminosities with 
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